Determining final probabilities directly from the initial state 
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In quantum scattering problems it is not possible to relate unambiguously a particular feature of 
the final outcome with some specific section of the initial state. As it is shown here, this drawback 
can be overcome by conveniently combining the divergence or Gauss-Ostrogradsky theorem with 
Bohmian mechanics. This renders a general approach which enables such a connection and allows 
us to determine the value of final partial or restricted probabilities directly from the corresponding 
localized section of the initial state. As an illustration, this approach is applied to two prototypical 
scattering phenomena: tunneling and grating diffraction. 

PACS numbers: 03.65.-w, 03.65.Ca, 03.65.Nk, 03.75.-b, 03.65.Xp 



I. INTRODUCTION 

Unlike classical mechanics, quantum mechanics poses 
an important problem when dealing with scattering prob- 
lems: although desirable, it is not possible to establish 
an unambiguous relationship between a particular fea- 
ture of the final outcome and some specific part of the 
initial state. Within standard quantum mechanics we 
find robust theoretical frameworks to tackle scattering 
problems [l| as well as efficient numerical algorithms Q 
when the former are not enough. However, they forbid 
the possibility to only compute the probabilities related 
to such features (i.e., partial or restricted probabilities) 
from a certain set of initial conditions, as it can be done 
in (statistical) classical mechanics by means of an appro- 
priate sampling. 

In this work we propose a way to attack this kind 
of problems, which combines the divergence or Gauss- 
Ostrogradsky theorem with Bohmian mechanics 0, Hj]. 
As it is well known, in standard quantum mechanics the 
divergence theorem allows to establish a connection be- 
tween the variation of the probability density p confined 
within a certain (configuration space) region of volume 
V with its flux across the boundary surface 5* of such a 
region, which is accounted for by the probability current 
density J. Formally, this is expressed in terms of the 
well-known continuity equation 



dp 
dt 



VJ = 0. 



(1) 



On the other hand, the possibility to define the analog 
of a classical trajectory in Bohmian mechanics, namely a 
Bohmian or quantum trajectory, is also of much help, for 
it allows us to associate ([1]) with ensemble-like descrip- 
tions. Thus, consider we want to compute the probability 
inside some subregion of the configuration space, E, that 
encloses one particular scattering feature. For example, 
this region may represent the boundaries embracing a 
particular diffraction peak in scattering problems or the 
transmission probability in tunneling problems. Usually, 
in a typical time-dependent quantum-mechanical simula- 
tion , one lets the process to propagate in time until a 
certain final configuration (state) is reached. In general, 



it is common to choose asymptotic, stationary configu- 
rations, which are closer to analytical time-independent 
derivations but this requirement is not always neces- 
sary. Then, one computes the partial fraction of the total 
probability that has gone into E, namely the restricted 
probability 0-0]) which is given by 



P E (oo) = lim Vs(t) 



lim 

t— >oo 



p(r,t)dr. 



(2) 



This quantity is precisely directly related to the problem 
posed above. 

As it is shown that, making use of Bohmian mechan- 
ics, if the set of quantum trajectories lying down on the 
boundary of E at t is known, restricted probabilities like 
© can be directly evaluated from the initial state at to. 
In virtue of the causal (Bohmian) connection between 
two configuration space points at two different times, 
the final positions of the trajectories forming the bound- 
aries of E are unambiguously associated with a certain 
set of initial conditions that constitute the boundaries 
of a region So at to- This implies that any trajectory 
starting inside So will end up inside E at t. Therefore, 
the restricted probability enclosed by the set of trajec- 
tories forming the boundaries of So will remain constant 
along time, allowing us to compute Ps(t) directly from 
the initial state. Indeed, it is also possible to obtain a 
one-to-one relationship between a specific part of the ini- 
tial state and a particular feature of the outgoing state 
that we might be interested in; at each time, one has 
full knowledge of where the probability confined within a 
particular So goes. It is worth noticing that, within the 
standard formulation of quantum mechanics, there is no 
way to determine this key information, which can be used 
in a rather practical way to conduct dynamical analyses 
of isolate scattering features in matter wave experiments 
Q or quantum control scenarios Q, for example. The 
generality of this result is proven for two prototypical 
processes found in scattering problems: tunneling and 
grating diffraction. 

This work is organized as follows. In Sec. [Til the gen- 
eral theoretical framework that justifies our approach is 
introduced. In Sec. IIII1 some applications will be con- 
sidered. First, as a simple analytical illustration of the 
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assertions discussed in Sec. [H] the case of a Gaussian 
slit diffraction will be exposed. Then, some numerical 
simulations for tunneling and grating diffraction will be 
presented. Finally, the main conclusions extracted from 
this work are summarized in Sec. IIVI 



II. THEORETICAL FRAMEWORK 

Let P{t) be the probability describing a certain prop- 
erty of a quantum system of mass m inside a volume V of 
its associated configuration space at a given time t, i.e., 



Pit) 



p(r,t)dr, 



(3) 



where p(r,t) = \^(r,t)\ 2 and the system wave function 
\1/ satisfy the corresponding Schrodinger equation. The 
variation with time of P(t) within V is given by 



dP(t) 
dt 



dp 
dt 



dr, 



which can also be expressed as 



dPjt) 
dt 



I (V • J) dr = - / J 

Jv Js 



dS, 



(4) 



(5) 



according to the continuity equation ([T]) , where dS in the 
second equality is a vector normal to a surface element 
dS on S and pointing outwards. From the application of 
the divergence theorem in this second equality, it follows 
that the losses or gains of P(t) inside V are described 
by the probability flow through S, which is quantified 
by the quantum current density J. Physically, this is a 
very clear and elegant result to describe the conservation 
of probability. However, from a more practical point of 
view, it presents a serious drawback. In order to under- 
stand it, consider a classical system, for which the same 
conservation equation holds (with p describing a density 
distribution). Let Ei and E2 be two regions with surfaces 
Si and 52, respectively. Both regions are assumed to be 
embedded in the system phase space at times t\ and t2, 
respectively, with t-i > t\. Moreover, each point on Si 
is causally connected with one point on S\ by means of 
a (classical) trajectory (solution of Hamilton's equations 
of motion). Then, any swarm of trajectories confined 
within Ei at t\ will end up inside E2 at £2, this being the 
trait of a conservative Liouvillian dynamics or, in other 
words, of the conservation equation for a classical sys- 
tem. In standard quantum mechanics, on the contrary, 
there is no way at all to define such a mapping transfor- 
mation, for there is no rule that allows us to connect the 
two regions Ei and E2. Even if one defines tubes along 
which the probability flows, there is nothing preventing 
it to scape from or to enter into such tubes, in general, 
because it is not possible to establish a causal connec- 
tion that describes the generation of these tubes when 
Ei evolves in time to E2. 



In order to establish such a connection one has to con- 
sider Bohmian mechanics, where the flow of quantum 
probabilities can be monitored by means of quantum tra- 
jectories (this flow, indeed, will define the probability 
tubes that we need). Bohmian or quantum trajectories 
are solutions of the equation of motion 



VS 

m 



J 

P 



(6) 



This equation is identical to the classical Jacobi law 
of motion [10], although here S refers to the phase of 
the wave function when it is written in polar form, i.e., 
^(r,^) = p 1 l 2 (v 1 t)e lSlyr ' t ^ h , instead of to a classical ac- 
tion. Nonetheless, it can also be defined if a local veloc- 
ity field is introduced according to the relation J = pv 
[second equality in ((B))], which describes how the quan- 
tum probability density is transported through the con- 
figuration space in the form of the quantum probability 
current density by means of a drift (local) velocity field 
v. Thus, consider Eo is a certain region of the system 
configuration space covered by swarm of initial condi- 
tions distributed randomly according to po( r ) = |^o( r )| 2 ; 
as it is also commonly found when dealing with statis- 
tical classical mechanics (i.e., just like in any standard 
classical Monte-Carlo sampling). The corresponding tra- 
jectories will then reach another region E at a time t 
causally connected with Eo, such that the positions of 
the trajectories follow a random distribution according to 
p(r,t) = |\I/(r, t) | 2 . Because of the non-crossing property 
of quantum fluxes (or, equivalently, quantum trajecto- 
ries) [Til ] . similarly to classical mechanics no (quantum) 
trajectories will scape from or enter into the probability 
tube defined by the evolution of Eo along time. There- 
fore, the probability inside, proportional to the number 
of trajectories confined within the boundaries of E, will 
remain constant along time. From now on, those trajec- 
tories starting just on points belonging to the surface of 
Eo will be called separatrix trajectories, for they separate 
the probability inside the tube from the rest in virtue of 
the non-crossing flux property. 

The previous argument can be reverted, this implying 
that given the end points of the separatrix trajectories 
associated with E, one can obtain the boundaries of Eo 
and, therefore, determine the probability inside E with- 
out carrying out the full simulation, but only knowing the 
associated section from the initial wave function (i.e., the 
portion of the initial wave function lying down within the 
boundaries of E ). In order to proof this statement, con- 
sider for simplicity (but without loss of generality) a one- 
dimensional problem, where E is defined as the region de- 
limited by two given points xf and x?? asymptotically. 
By "asymptotic" it is meant the Fraunhofer domain |12l ] , 
where the probability density remains stationary, i.e., its 
relative shape does not change with time, but remains 
re-scalable. As shown in Ref. [lj], these points can be 
associated with the positions of two boundary (separa- 
trix) quantum trajectories (which, within the Fraunhofer 
regime, are also re-scalable, following the laws of uniform 
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motion). Because these trajectories can be propagated 
backwards in time according to Eq. ©, Eq. (J5|) can be 
reexpressed as 



p(x°°)dx c 



x 2 (t) 



= lim / p[x{t)]dx{t) = lim T>s{t), (7) 

where p(x°°)dx 00 is the probability to find the system 
confined between x°° and x°° + dx°°; p[x(t)]dx(t) is the 
probability to find the system confined between x(t) 
and x(t) + dx(t); and the time-dependence denotes a 
causal connection, e.g., x°° — \\m t ^oa x{t). The time- 
dependence of the probability density is purposely spec- 
ified as p[a;(i)] instead of p(x,t) in order to stress the 
fact that, at any time t, it is evaluated along quantum 
trajectories rather than on arbitrary points of the config- 
uration space (this also emphasizes the fact that quan- 
tum trajectories are paths along which probability flows 
in configuration space). 

Since Bohmian trajectories can be traced back in time, 
from t to their initial condition, the limit in the second 
equality in Eq. (J7]) can be removed, so that V\z{oo) can 
be recast as an integral over the system initial state, 



Pe(oo) 



x 2 (0) 



Ei(O) 



p[x(0)]dx(0) =Vs(0), 



(8) 



which is the same to say that it can be obtained just 
from a sampling over initial conditions (for the quan- 
tum trajectories). This is a result of general validity, 
which just states that the probability in a certain region 
of the configuration space can be transported to another 
one unambiguously, i.e., following causally (by means of 
Bohmian mechanics) the evolution of the boundaries of 
such an initial region. As it can be inferred, this result 
goes beyond standard quantum mechanics, for although 
it is possible to define tubes along which the probability 
flows and surfaces traversed by such flows, as mentioned 
above there is no rule accounting for the evolution of 
probability tubes, thus avoiding to establish any univo- 
cal relationship between regions £ defined at different 
times. As a consequence, we find that, in principle, the 
knowledge of both the initial state and some particular 
initial conditions of interest, .ti(0) and £2(0), allows us 
to determine quantities like transmission probabilities or 
diffraction intensities (or an estimate of them without 
performing the whole dynamical calculation provided we 
have a fair guess [H|]). 

Furthermore, it is also worth stressing the strong con- 
nection between Eq. ([5]) and the so-called Born rule [3] . 
Indeed, the combination of the quantum continuity equa- 
tion and Bohmian mechanics makes Eq. to incorpo- 
rate implicitly a kind of time-dependent Born rule. In 
this regard, it is interesting to note that, instead of the 
two limiting times t — > 00 and t — 0, one can also con- 
sider two arbitrary times, t and t' (we will assume t' > t). 



Then, from Eq. (JSJ, it follows that 



a(t') px 2 (t) 

p[x(t')]dx(t') = / p[x{t)]dx(t). (9) 

xi{f) Jxi(t) 



On the other hand, because of the causal connection be- 
tween x(t) and x(t') in Bohmian mechanics, one can also 
define a Jacobian 



J[x(t)} 



dxjf) 
dx{t) 



(10) 



which describes the mapping transformation in configu- 
ration space from x(t) at a time t to x(t') at a time t' , as 
it is well known from classical mechanics [HI, HH . Tak- 
ing thus into account the equality (|9|) and the connection 
between dx(t) and dx(t') enabled by the Jacobian, 



dx(t') = \J[x(t)]\dx(t), 



(11) 



one finds that the probability density evaluated along a 
Bohmian trajectory at a given time t' is related through 
the inverse Jacobian transformation with its value at an 
earlier time t, i.e., 



p[x(t>)] = \J[x(t>)}\p[x(t)}, 



(12) 



with |v7[z(t')]| = iJ'^Wir 1 . In other words, this result 
states that Born's rule preserves along time whenever the 
evolution of the probability p[x(t)] is monitored between 
two Bohmian trajectories x{t) and x'(t) = x(t) + dx(t), 



i.e. 



p[x(t')]dx(t') = p[x(t)]dx(t), 
as in classical mechanics [TH, [lj| . 



(13) 



III. APPLICATIONS 

A. Gaussian slit diffraction 

The free evolution of a quantum system of mass m 
after undergoing diffraction by a Gaussian slit constitutes 
a simple example to illustrate the assertions from Sec. ITT] 
Thus, consider the wave function describing the system 
just after the diffraction process is given by 



(14) 



where Aq = (27T(Tq) _1 / 4 is the normalization constant, 
xq and po are respectively the (initial) position and mo- 
mentum of its centroid [i.e., (x) = xq and (p) = po\, and 
do its initial spreading. The time-evolution of this wave 
packet is described by 

\£( x f\ — A t e -(x-x cl ) 2 /i5- t o-o+ipci(x-x cl )/h+iE cl t/h j-^g-j 

with A t = (27ro- t 2 ) -1 / 4 , a t = a (l + iht/2ma$) and 
E c i = Po/2to. The shape of this wave packet does not 
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change with time, but only its width undergoes a mono- 
tonic increase, as 



0* = Ft 




(16) 



while its centroid evolves according to a classical uniform 
motion, i.e., x c \ = x + v t and p c \ — po- The associated 
probability density thus looks the same at any time, 



p{x,t) = 



1 



,-(x-x*) 2 /2a* 



(17) 



The Bohmian trajectories associated with (p~5]) are readily 
obtained after integrating the corresponding equation of 
motion, to yield 



x(t) = x c i + — [x(0) - x ] 

CO 



(18) 



From this expression, the associated Jacobian reads as 

dx(t) O t 



J[x(0)} 



dx(0) 



oo 



which leads to 



dx{t) 
p(x(t)) 



\J[x(0)]\dx(0), 
\J[x(t)}\p(x(0)), 



(19) 



(20) 
(21) 



with \J[x(t)}\ = |JX0)]| _1 = <T /<T t ruling the back- 
wards mapping transformation, from t to to — 0. If these 
expressions are now substituted into Eq. ([7]). we readily 
find Eq. JgJ or, equivalently, Eq. (fig]). 
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FIG. 1: (a) Initial (black line) and asymptotic, final (red 
line) probability densities representing the tunneling of a wave 
packet with an arbitrary shape through a quasi-squared bar- 
rier (blue shadowed region), (b) Time-dependence of the 
different restricted probabilities: transmission (green dashed 
fine), reflection (blue dash-dotted line) and intra-barrier reso- 
nance (red dotted line) probabilities. The transmission prob- 
abilities directly computed from the initial state accounting 
only the associated quantum trajectories is displayed with 
black solid line, (c) Quantum trajectories illustrating the tun- 
neling dynamics. The separatrix trajectory is denoted with 
the red thicker line, (d) Separation in two sections of the ini- 
tial probability density according to the separatrix trajectory. 
Only the green shadowed region will contribute to the trans- 
mittance, while the blue shadowed section will be scattered 
backwards, contributing to the reflectance. 



B. Tunneling 

Usually, within the energy domain, the calculation of 
properties related to tunneling through barriers (trans- 
mittances, reflectances or intra-barrier resonances) is car- 
ried out by means of stationary wave approaches, where 
the system wave function and its first derivative have to 
satisfy some matching and/or boundary conditions at the 
barrier edges [H, [I3|- The approach presented above, on 
the contrary, enables the calculation of such properties 
in a more straightforward fashion, directly from the ini- 
tial state once the separatrix trajectories are determined 
flijj . Thus, we have considered the scattering with a bar- 
rier of a wave packet consisting of a coherent sum of three 
Gaussian wave packets, 

3 

*o(a0 =^2ciipi{x), (22) 

i=l 

with each ipi having the same form as (|14l) . In 
particular, in the example shown here the pa- 
rameters chosen are: (ci, 02,03) = (1.0,0.75,0.5), 

(xq,i,x ,2, 2:0,3) = (-10,-12,-9), (po,i,Po,2,Po,3) = 



(10,20,15), (<7 ,i, 0-0,2, 0-0,3) = (0.2,1.6,0.8), and m = 
H = 1; after introducing these values, the wave function 
^0 was properly renormalized before starting the sim- 
ulation. Regarding the barrier, it was a quasi-squared 
barrier consisting of a sum of two hyperbolic tangents, 

V(x) = — {tanh [a(x — X-)\ — tanh [a(x — x+)]} , 

(23) 

with Vq — 150, a = 10, and x± — ±2. As shown below, 
this example will allow us to illustrate the generality of 
Eq. ((5J with respect to the election of the initial state 
(initial condition) or the type of barrier (problem) . 

The initial and final probability densities are displayed 
in Fig.QJa) (with black and red solid lines, respectively); 
the barrier has also been plotted (blue shadowed region). 
Regarding the final state, the system has been let to 
evolve until the probability in the region covered by the 
barrier (intra-barrier resonance) was negligible, this be- 
ing considered the asymptotic time. This can be better 
appreciated in Fig. Hlb), where the transmission (green 
dashed line), reflection (blue dash-dotted line) and intra- 
barrier resonance (red dotted line) probabilities are mon- 



5 



itored along time. These three restricted probabilities are 
computed assuming the corresponding regions are as fol- 
lows: Straus is the region beyond the right edge of the 
barrier, E rcs is the region bound between the two edges 
of the barrier, and E re fl is the region to the left edge of the 
barrier. As it can be seen, after t rj 1.15, Vs IBB « and 
Ps trans reaches its maximum, asymptotic value; "Px; rcfl , on 
the contrary, reaches its minimum earlier, at t » 0.75. 

Taken that into account, we decide to split the asymp- 
totic state between transmission and reflection, with the 
region So being being associated with the former. Ini- 
tially, the upper boundary of this region can be any tra- 
jectory to the right border of the initial probability den- 
sity, such that the associated probability density is neg- 
ligible, i.e., p(xo) w 0. For the lower boundary, a search 
has to be done [13]], so that it is ensured that the cho- 
sen trajectory is the last one in passing the barrier with 
a negligible value of the associated probability density. 
This can be determined by means of a series of sampling 
runs or just fixing the asymptotic value of the trajec- 
tory and running backwards the dynamics until t = 0. 
In Fig. HJc), a sampling set of quantum trajectories is 
shown. These trajectories have been evenly distributed 
along the whole extension covered by the initial proba- 
bility density to better illustrate the dynamics. In this 
case, the red thicker line denotes the separatrix trajec- 
tory which separates the initial swarm into two groups: 
those that will surmount the barrier and those that will 
bounce backwards. This means that, in this case, the 
region E (the time-evolved of Eq) is specified only by 
one separatrix trajectory. When this is translated into 
the initial probability density, we find which part of it 
will contribute to transmission and which one not [see 
Fig. Did)]. Thus, if we compute the integral over the 
green shadowed region displayed in Fig. QJd) and follow 
it along time taking into account where the separatrix 
trajectory is at each time, we obtain the black solid line 
plotted in Fig.QTc). This curve presents a constant value 
at any time, which coincides with the asymptotic value of 
the transmission probability, thus proving our assertion 
that final probabilities can be computed directly from 
the initial state and, as it will be shown with the next 
example, this can only be done by means of Bohmian 
trajectories. 



C. Grating diffraction 

Quantum trajectories have been used in the past to 
determine the boundaries in the initial probability den- 
sity that separate contributions associated with different 
specific features in both [l2|] and slit diffraction [l^]. By 
means of the approach introduced in Sec. [Ill one can take 
these studies further away, to a more quantitative and ap- 
plied level. Thus, consider one would like to determine 
an observable in any of these scattering problems, such 
as the so-called peak-area intensity, for example, which is 
the integral of the probability density covered by a partic- 



ular diffraction peak. Physically, this quantity provides 
us with the relative amount of diffracted particles lying 
within the width covered by the corresponding diffraction 
peak. According to the tools developed so far, in partic- 
ular Eq. (|l?p. this quantity can be obtained directly from 
the initial wave function once the two trajectories delim- 
iting the corresponding peak are determined, just as in 
the example of tunneling described above. To illustrate 
this problem, for simplicity (but without loss of gener- 
ality), we consider the case of multiple-slit diffraction, 
which is isomorphic to that of scattering by a periodic 
reflection grating [3]. Thus, a five-slit array of Gaus- 
sian slits is assumed, so that the outgoing wave can be 
represented by a coherent superposition of five identical 
Gaussian wave packets, like (fl4l . 

5 

*„(*) =X>(*), (24) 

i=l 

where Xq = —4 + 2(i — 1), p$ = po = (zero transverse 

momentum), = (To = 0.2, and m = h = 1; as before, 
after introducing these values, the wave function was 
properly normalized before conducting the simulation. 
Note that we have implicitly assumed that the propaga- 
tion along the direction perpendicular to the grating is 
too large compared to the evolution along the transversal 
one, this enabling this simplified one-dimensional treat- 
ment, where we only focus on the diffractive and inter- 
ference effects that take place along the transversal di- 
rection. Furthermore, now the wave packet will not be 
allowed to evolve asymptotically, until all the minima are 
well defined (i.e., when the probability density at such 
values vanishes) , but only until they can be well resolved 
(i.e., when the distance betweeen two maxima is such 
that they can be well distinguished), as it can be seen 
in Fig. [2a), for a propagation up to t = 10. This may 
be the case when the detector is allocated relatively close 
to the diffractive object, so that the system has not had 
time enough to reach the Fraunhofer regime, where the 
diffraction pattern becomes stationary (its general profile 
does not depend on time). In such a case, we can ask our- 
selves questions like which parts of the diffracted beam 
contribute to some particular diffraction peak, or what is 
the contribution of a particular slit to the final intensity 
collected (or, equivalently, to one particular diffraction 
feature). 

Consider the principal maxima associated with the 
diffraction orders n = and n = +1 in Fig. (UJa). We de- 
fine Eo = {xq~,2;J} and E+i = {#Xi,2fti} as the regions 
bound by the adjacent minima to each one of these max- 
ima, respectively, at t — 10 [see color shadowed regions 
in Fig.(2Ja)]. The restricted probabilities associated with 
these two regions give us the value of the correspond- 
ing peak-area intensities. However, unlike the tunneling 
problem, the calculation of the time-evolution of these 
probabilities from a standard quantum point of view is 
here a subtler issue, for there is not a specific space re- 
gion that can be uniquely defined at any time. Only once 
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FIG. 2: (a) Interference pattern obtained from the superposi- 
tion of five Gaussian wave packets at t = 10. Principal max- 
ima are labeled according to their diffraction order n and the 
color shadowed regions denote for which the peak-area inten- 
sities are computed, (b) Peak-area intensities associated with 
the diffraction orders n — (blue line) and n = +1 (red line). 
The peak-area intensities directly computed from the initial 
state and using the separatrix trajectories at t — 1 are dis- 
played with black solid line in both cases. Black dashed lines 
indicate the asymptotic value of the peak-area intensities, (c) 
Set of quantum trajectories leaving each slit. The straight 
color lines delimit the tubes of probability considered to com- 
pute the peak-area intensities of the corresponding colors in 
panel (b). (d) Separation of the two sections in the initial 
probability density according to the separatrix trajectories at 
t — 10. The blue shadowed section refer to n = 0, while the 
red one to n = +1. (e) Contributions from slits 3, 4 and 5 
[see panel (d)] to the different diffraction orders along time, 
(f) Set of quantum trajectories around the separatrix trajec- 
tory at t — 10 (purple thicker line). The right-edge boundary 
of the n — probability tube is denoted with the straight 
blue thicker line, while a separatrix trajectory at a later time 
(t ~ 20) is denoted with green thicker line. 



the maxima can be resolved (i.e., we observed a minimum 
between two adjacent maxima), we can consider a bound- 
ary, but not earlier, when the diffraction pattern is still 
in the Fresnel regime and there is no clue on which part 
of the pattern can be associated with a particular fea- 
ture in the far field or Fraunhofer regime. Thus, the only 
thing we can do is to define a series of angular sectors or 



probability tubes, which start at the center of the grating 
and prolong (in time, because this is a one-dimensional 
problem) up to the position of each adjacent, vanish- 
ing minimum of a given diffraction (principal) intensity 
peak. That is, if such adjacent minima appear respec- 
tively at positions x^ at a time t for the maximum of 
order n, the time-evolution of the boundaries of the cor- 
responding tube will be given by cc„ (t) = (x^ /t)t, i.e., at 
t the peak-intensity will be computed between x~ (t) and 
x+(t). Note that this hypothesis is based on assuming 
the wave function is already well inside the Fraunhofer 
regime, where the minima are well defined (i.e., they are 
zero) and therefore the probability inside the tubes will 
remain constant because has already reached the station- 
arity (the relative shape does not change with time; time- 
dependent scattering approaches are based precisely on 
this idea [![). This asymptotic trend can be observed in 
Fig. [D^b) , where the peak-area intensities for the diffrac- 
tion orders n = (blue dashed line) and n — +1 (red 
dotted line) are plotted against time. As it can be seen, 
these quantities increase with time until they start to ap- 
proach an asymptotic value, denoted by the black dashed 
lines. These lines denote the peak- area intensities ob- 
tained when the minima are zero and, in virtue of the 
non-crossing property of quantum fluxes [ll[ , no transfer 
of probability is possible between adjacent maxima. In 
particular, at t = 20 the deviations with respect to the 
asymptotic values are about 1.72% for n — and 1.68% 
for n = +1, which are essentially negligible. 

As in the case of tunneling, if we now consider the tra- 
jectories that reach the minima adjacent to the peaks 
n = and n = +1, we can unambiguously delimit 
the regions Eq and E^ (the subscript "B" denotes the 
Bohmian simulation) that allow us to compute from the 
beginning the corresponding restricted probabilities. In 
Fig. [2c) we observe a set of trajectories illustrating the 
dynamics of the diffraction process; to compare with, the 
thicker blue and red boundaries of the tubes binding Eo 
and E + i, respectively, have also been displayed. By fol- 
lowing backwards the separatrix boundaries for Eq and 
E?!, we find the initial positions that allow us to know 
which portions of the initial wave function will contribute 
to each one of the two features considered, as shown in 
Fig. Hfd). In this figure we note that the central maxi- 
mum is formed by the accumulation of particles that left 
the central slit and almost a half of the two neighbor- 
ing slits. On the other hand, the maximum n = +1 is 
built up by a small portion of particles leaving the slit 
centered at x = 2 and most of those coming from the 
slit centered at x = 4. The contributions arising from 
each slit and each particular region of each slit remain 
constant along time, as shown in Fig. HJe). Again, this is 
due to the non-crossing property, which avoids trajecto- 
ries inside a certain tube of probability to abandon it or 
to allow others from outside to enter in it, thus keeping 
the number of them constant along time. Taking this into 
account, we find in Fig.^b) that the restricted probabili- 
ties inside the boundary trajectories displayed in Fig.[5Jc) 
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remain constant along time and, at t = 10, they corre- 
spond exactly with those that one would expect from the 
peak-area calculations (blue and red curves); the cross- 
ing between the standard and the Bohmian calculations 
is indicated by the downwards pointing arrows. As it 
can also be noted, beyond t = 10 these contributions are 
smaller than the peak- area ones, which is due to the fact 
that we do not have a clear distinction in terms of trajec- 
tory conditions between those trajectories that will go to 
the central maxima and those remaining in the adjacent 
secondary maxima. 

As it can be seen in Fig.^f), which is an enlargement 
of panel (c), as time proceeds there are still trajectories 
getting into the sector associated with the maximum n — 
0. In this panel we observe how the right-side separatrix 
trajectory for t ~ 20 (green thicker line) appears a bit 
further away than that for t = 10 (purple thicker line). 
Although the difference between the initial condition xq 
for these two trajectories is only about 0.030 (i.e., it is 
negligible in practice), as time proceeds the difference 
is becomes rather significant. This implies that all the 
trajectories between the purple and the green ones at 
t = 10 will not be considered in the calculation of V, 
but only gradually at subsequent times. Of course, if we 
compute V taking into account the purple line all the 
way (until t = 20), all those trajectories between it and 
the boundary of the probability tube [blue thicker line in 
Fig. |U[f)] will not be taken into account and, therefore, 
the value of V will be slightly lower than the asymptotic 
value V n [oo). For example, here we find that Vq and V+\ 
are about 6.92% and 6.70% lower than the corresponding 
asymptotic values, which means that a 6.92% and an 
6.70% of trajectories are still lacking in the calculation 
of the corresponding peak areas. 



IV. CONCLUSIONS 

Here we have presented a method that allows us to re- 
late unambiguously final probabilities with different sec- 
tions of the initial state, which is closely connected to the 
standard methods that can be found in statistical clas- 
sical mechanics, e.g., in atmospheric modeling (l9l. [20|. 
The possibility of such a connection arises from the com- 



bination of the divergence or Gauss-Ostrogradsky theo- 
rem with Bohmian mechanics, which allows us to go a 
step beyond the standard formulation of quantum me- 
chanics and devise new tools to study and to understand 
the physics underlying of quantum systems. In particu- 
lar, even if not observable from a experimental point of 
view, the fact that one can properly define (in terms of 
separatrix trajectories) regions which are uniquely trans- 
ported throughout configuration space without loss or 
gain of quantum flux results of much interest to study 
the dynamical role of quantum phase in time-resolved 
experiments. Note that this methodological prescription 
allows to monitor the detailed evolution of a set of par- 
ticular initial conditions and, therefore, to determine the 
final outcome (probability) directly from the initial state. 
Actually, even if one does not know where such a partic- 
ular set or part of the initial wave function will evolve, 
it is sure that the probability confined within its bound- 
aries (separatrix trajectories) will not be mixed up with 
contributions from other parts of the same initial wave 
function. Within the standard scenario, however, there is 
no certainty about this and, therefore, one has to appeal 
to rather arbitrary methods to determine final restricted 
probabilities. 

The approach posed here was aimed at computing 
probabilities without any need for solving the full dy- 
namics of the process, but only with some knowledge pro- 
vided by the quantum trajectories. In the literature it is 
also possible to find methods that operate the other way 
around, i.e., from probability densities they try to infer 
the corresponding quantum trajectories without solving 
the associated equations of motion. This is the case, for 
example, of the Bohmian Monte Carlo sampling method 
21 1, based on the idea of quantile motion [22||, or the 
kinematic approach [23j , based on the Voronoi's tessella- 
tions method (2~3 |. 
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